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This  report  describes  some  simple  techniques  for 
smoothly  filling  in  gaps  in  object  contours.  The  first 
technique  considered  was  recently  proposed  by  Ullman  [I] ; 
it  constructs  the  completion  of  the  contour  using  two  arcs 
of  circles  that  are  tangent  to  the  gap  ends  and  to  each 
other,  and  that  have  minimum  total  curvature.  An  analysis 
of  this  technique  is  presented,  and  exeunples  of  its  use 
are  given.  A second  technique  uses  cubic  polynomial  com- 
pletions; when  suitably  constrained,  this  technique  yields 
very  reasoneUale  completions. 
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This  report  describes  some  simple  techniques  for  filling 
in  gaps  in  object  contours.  Such  gaps  occur  in  many  classes 


of  images  as  a result  of  occlusion,  shadows,  etc.  Many 
levels  of  knowledge  may  be  broght  to  bear  upon  the  problem  of 
filling  in  a missing  segment  of  an  object  boundary.  The  tech- 
niques considered  here  utilize  only  local  information. 
Specifically,  we  are  given  two  points  through  which  the  con- 
tour must  pass,  and  the  directions  in  which  the  contour  is 
going  at  those  points. 


1.  The  Ullroan  Method 

The  first  technique  we  will  consider  was  proposed  by 
S.  Ullman  (!];  the  following  description  is  quoted  from 
Ullman's  paper: 

Subjective  contours  are  perceived 
when  the  visual  system  fills  in  the  gap 
between  distinct  edges.  Other  examples 
of  filling-in  processes  are  the  perceived 
trajectory  in  optimal  Beta  motion  and  the 
continuation  of  the  visual  field  across 
retinal  scotomata.  The  study  of  these 
filling-in  processes  has  concentrated  to 
date  on  the  triggering  conditions.  For 
instance,  under  what  conditions  will  a 
complete  trajectory  be  perceived  between 
successively  presented  stimuli,  or  what 
conditions  will  enhance  the  generation  of 
subjective  contours  between  distinct 
edges.  Another  important  but  hitherto 
neglected  problem  posed  by  the  filling-in 
phenomena  concerns  the  shape  of  the  filled- 
in  contours  and  trajectories.  Since  sub- 
jective contours  are  synthesized  by  the 
visual  system  the  excunination  of  their 
I shape  might  provide  clues  about  the  nature 

; of  the  mechanism  that  generates  them.  In 

I this  paper  I shall  first  excunine  the  shape 

I of  subjective  contours  and  then  suggest  a 

\ ^ network  capable  of  generating  these  shapes. 

I Of  the  infinitely  many  loci  passing 

[ through  two  given  edges,  which  is  the  one 

j chosen  by  the  visual  system?  Let  me  start 

I the  quest  for  an  answer  by  stating  three 

j observations  and  an  hypothesis  which  will 

i 


i 


1 


serve  as  guidelines  for  narrowing  down  the 
range  of  possible  shapes  to  a unique  one. 

The  resulting  curves  will  then  be  compared 
» h actually  perceived  filled-in  contours. 

I.te  four  guidelines  are: 

1.  Isotropy;  As  far  as  the  filled- 
in  contours  are  concerned  the  visual  field 
seems  approximately  homogeneous  and  isotropic: 
The  filled-in  contours  produced  by  a given 
figure  do  not  change  in  shape  when  the  figure 
is  translated  or  rotated. 

2.  Smoothness:  except  for  some  special 
cases  of  filled-in  corners,  the  generated 
curves  are  smooth,  that  is,  dif ferentiadile 

at  least  once. 

3.  Minimum  curvature:  This  guideline 
is  inspired  by  the  resemblance  of  the  filled- 
in  contours  to  a thin  doubly  cantilevered  beeun 
or,  alternatively,  to  the  curve  known  in 
approximation  theory  as  a cubic  spline.  Both 
curves  are,  in  some  sense,  loci  of  minimum 
curvature.  The  shape  of  a subjective  contour 
passing  through  two  given  edges  whose  orienta- 
tion difference  is  much  less  than  7r/2  closely 
resembles  the  cubic  spline  passing  through  this 
boundary  edges.  At  this  point  I only  wish  to 
state  this  resemblance  informally:  the  minimum 
curvature  property  will  receive  a formal  treat- 
ment in  the  sequel. 

4.  The  locality  hypothesis;  The  operation 
by  which  a boundary  edge  is  extended  to  generate 
a subjective  contour  is  assumed  to  be  local  in 
nature.  That  is,  it  depends  only  on  the  end 
part  of  the  given  edge  to  be  extended  and  not 

on  the  shape  of  the  entire  edge. 


The  four  guidelines  are  sufficient  to  resolve 
the  shapes  of  the  filled-in  contours. 

Let  the  set  of  possible  extensions  (of 
a contour]  A be  LI,  12,..., Lk.  The  set  of 
curves  Ll,...,L)c  possesses  the  following  pro- 
perty. If  you  take  any  curve  Li  and  cut  it 
at  scxne  point  P,  then  the  shape  of  P-Li,  the 
continuation  of  Li  from  P onwards,  is  con- 
gruent to  one  of  the  curves  Ll,...,Lk.  One 
can  meet  this  requirement  by  taking  all  the 
Li's  to  be  arcs  of  some  circles  (including 
the  limit  case  of  a straight  line) , since 
any  part  of  a circle  is  congruent  to  any 
other  part  of  equal  length.  Under  certain 
conditions,  the  circles  solution  is  the  only 
set  of  carves  obeying  the  above  requirement; 
however,  I shall  not  diverge  here  to  examine 
these  conditions. 

Since  the  produced  curves  must  all  be 
tangent  to  the  edge  A to  meet  the  smoothness 
requirement,  the  prediction  is  that  the 
filled-in  contour  is  composed  of  the  arcs  of 
two  circles,  one  temgent  to  one  edge,  the 
other  tangent  to  the  second  edge.  An 
additional  restriction  is  imposed  by  the 
smoothness  guideline,  which  requires  that  the 
two  arcs  share  a ccxnmon  tangent  at  their 
meeting  point.  There  is  still,  however,  an 
infinite  niutiber  of  arc  pairs  satisfying  all 
the  above  conditions.  Following  the  minimxim 
curvature  hypothesis  let  us  pick,  out  of  all 
the  admissible  pairs,  the  pair  which  minimizes 
the  total  curvature.  The  total  curvature  of  a 
curve  is  defined  here  as  the  integral  / (do/dl) 
over  the  curve,  where  a is  the  slope  of  the 


curve.  The  value  of  da/dl  at  a given  point 
P is  Icnoim  as  the  local  curvature  at  that 
point.  In  the  sequel,  hovrever,  the  word 
"curvature"  will  refer  to  total  curvature 
unless  otherwise  stated. 

Ullman  goes  on  to  describe  a locally  connected  network 
which  computes  these  contours.  In  the  absence  of  such  a 
network,  a conventional  approach  must  be  employed  to  compute 
the  contour.  We  will  now  describe  such  an  approach. 


Implementation  of  the  U liman  Method 


In  Figure  la  we  have  a circular  arc  emanating  from  point 
P in  direction  and  terminating  at  point  Q in  direction 
The  ^mgle  a is  measured  from  the  initial  direction  vector 
to  a vector  pointing  from  the  initial  point  of  the  arc  to  the 
terminal  point.  Note  that  the  total  angular  change  (f  d0)  is 
2a.  This  holds  true  for  o cf-n,  + it]  . We  run  into 
difficulty  if  we  use  a different  range  for  a.  For  example, 
in  Figure  lb  we  have  a * -ti/2,  and  the  arc  goes  through  an 
angular  change  of  2a  = -ir.  However,  if  we  instead  measured 
a as  3Tr/2,  and  constructed  a circular  arc  with  total  angular 
change  2a,  we  obtain  the  situation  shown  in  Figure  Ic.  The 
vector  from  the  initial  to  the  terminal  point  of  the  arc  now 
points  in  a direction  opposite  to  the  direction  defined  by  the 
angle  3ir/2.  Of  course,  it  is  obviously  foolish  to  construct 
completions  using  circular  arcs  with  angular  changes  exceeding 
2tt  in  magnitude.  Therefore  we  can  restrict  2a  to  (-2Tr,  +2ir)  , 
or  at  (-n,  +ir)  , for  the  arcs  used  in  our  contours. 

It  is  a matter  of  elementary  geometry  to  show  that  for 
Figure  la, 

f,  d©  2 4asina 

di  ^ D — • 

We  can  now  define  our  problem  with  the  aid  of  Figure 
2a.  The  contour  we  wish  to  compute  consists  of  two  circular 
arcs  which  join  at  point  J.  For  angles  cp^^,  ©j^,  and  a^^, 
associated  with  the  point  Pj^,  let  us  consider  counterclockwise 
rotation  as  positive,  while  for  cp2 » ®2'  “2'  will  consider 




clockwise  rotation  positive.  Although  this  may  seem  confusing 
at  first,  it  leads  to  symmetry  in  the  equations  developed  be- 
low. Thus  in  Figure  2a,  9^^,  (p^,  9^,  <p2  positive,  and 

and  02  would  be  negative.  All  angles  are  in  [-ir,  -fn]  . 

Let  us  now  consider  the  total  angular  change  in  follow- 
ing the  contour  from  to  Pj.  In  going  from  Pj^  to  point  J, 
the  angular  change  (considering  counterclockwise  as  positive) 
is  2a^,  as  discussed  earlier  in  connection  with  Figure  1.  As 
we  go  from  J to  P^,  the  angular  change  (again  considering 
counterclockwise  as  positive)  is  202 • At  first  this  may  seem 
puzzling,  since  02  is  measured  considering  clockwise  rotation 
as  positive.  Indeed,  the  counterclockwise  rotation  would  be 
-202  if  we  were  proceeding  from  P2  to  J;  however,  since  we  are 
moving  from  J to  P2,  the  counterclockwise  rotation  is  202- 
Thus  the  total  angular  change,  travelling  from  P^^  to  P2,  is 
2a^  + 202,  where  the  angular  change  is  measured  according  to 
the  conventions  used  at  P^.  The  total  angular  change  from  P^ 
to  P2  is  also  given  by  -^2  (again  measured  with  respect  to 
the  conventions  at  P]^)  • For  exaunple,  in  Figure  2b  the  total 
average  change  from  P^^  to  P^  is  -n.  Thus  we  have 
2oj^  + 202  =~('¥]^+'P2^  * course,  we  can  add  or  subtract 

multiples  of  2tt  to  the  total  angular  change  required  to  go 
from  Pj^  to  P2»  so  that  in  general 

20j^  + 202  = -(cpj^+92^ 


cpj^+cp2 

or,  defining  ^ — 2 — 


r? 

I. 

i' 

Therefore,  since  “ ‘J’l  "*■  ®2  " ***2  * **2'  have 

+ ©2  » ♦^yQ  + nn.  Since  ♦,  and  *2  ^ I "*»+*!»  clearly 
*AVG  fact,  can  only  equal  tn  if  both 

emd  *2  * -’’•  this  case,  illustrated  in  Figure  2c,  the  en- 
tire technique  breaks  down,  as  there  is  obviously  no  pair  of 
circular  arcs  which  will  join  Pj^  and  P2  as  required.  Hence, 

*AVG  situation  in  which  it  is  possible  to 

construct  the  contour.  Now  if  we  wish  to  form  a triangle 
Pj^  J P2  as  in  Figure  la,  then  |0j^  + ©2!  * n.  Thus: 

n can  equal  +1  only  if  ♦^yQ  * 0 
n can  equal  -1  only  if  ♦^yg  ^ 0 
|n|  can  never  exceed  1. 

To  summarize,  ©j^  and  ©2  obviously  determine  the  point  J at 
which  the  circular  arcs  join,  and  we  have  now  derived  a con- 
straint on  ©j^  and  ©2.  It  is  a constraint  which  has  a simple 
geometric  interpretation.  The  locus  of  points  J for  which 
©j^  + ©2  = ‘**avg  ^ ” restricted  as  given  above,  is  a 

circle.  For  example,  in  Figure  3 we  show  cases  where 

= Tf/Z,  Tf/ii  and  0.  The  last  of  these  cases  is  unusual  in 
that  we  can  only  form  degenerate  straight  line  triamgles. 

Although  it  is  possible  to  compute  contours  using  the  con- 
straint ®2.  ^ ®2  ” ^AVG  ^ with  non-zero  n,  we  can  get  into 
serious  difficulties  when  doing  so.  Note  that  we  are  attempt- 
ing to  find  a contour  which  fits  the  above  constraint  and 
has  minimum  curvature.  Now  consider  the  case  shown  in  Figure 
4a.  If  we  look  at  the  family  of  completions  defined  by 

I 


r 


kk. 


®1  ®2 


AVG 


0,  we  see  that  they  look  as  shown  in  Figure 


4b,  and  the  minimum  curvature  contour  obeying  the  constraint 
9]^  @2  = 0 is  as  shown  in  Figure  4c,  as  we  would  expect  from 

symmetry.  Now,  however,  consider  the  family  of  contours  obeying 


the  constraint  0j^  ®2  These  contours  look  as  shown  in 


Figure  4d.  The  curvature  of  a circular  arc  is 


2 

/(  ^ ) dZ  = ^ where  2a  is  the  angular  change  along  the  arc. 


d£  ' R 

and  R is  the  radius  of  curvature.  Note  that  in  the  family  of 
contours  illustrated  in  Figure  4d,  and  02  remain  fixed  at 
n/2,  while  the  radii  can  increase  without  bound.  Thus  there  is 
no  finite  solution  in  the  case  * 0,  or  «pj^  « ”***2*  This 

special  symmetric  case  with  »*  9^  is  not  the  only  problem. 
Consider  the  "almost  symmetric  case"  in  Figure  5a.  The  locus 
of  points  at  which  the  arcs  can  join  is  illustrated  in  Figure 
5b.  If  we  examine  the  feimily  of  contours  obeying 
0,  + 0-  = we  get  a minimum  solution  which  is  "almost 

symmetric"  (Figure  5c) . However,  there  exist  solutions  obeying 


0j^  + 02  = “■'f  which  have  lower  curvature  than  the  curve 


above  (assuming  close  enough  to  0) , as  in  Figure  5d. 

Although  in  this  case  ^ 0»  so  we  can  find  a finite 

minimum  curvature  contour  obeying  "*■  ®2  ~ ***AVG  “ find 

that  the  shape  of  the  contour  is  not  "almost  symmetric"  as  we 
would  intuitively  like  it  to  be.  In  fact,  the  closer  we 
approach  the  symmetric  case  = ~92'  larger  the  minimum 

contour  obeying  ®2  “ ^AVG  grows. 


-IT  grows. 


Because  of  this,  we  shall  from  now  on  use  + 9^  **  ^aVG* 

Remembering  again  that  8^^  and  are  angles  in  the  triangle 
?!  J P2'  it  is  clear  that  0j^  and  ©2  must  have  the  same  sign 
or  we  cannot  form  the  triangle.  Therefore,  ©j^  and  ©2  can  only 
take  on  values  between  0 and  There  is  one  further  point 

to  be  checked:  “ ®i  “ ‘J’l  * ®2“  ^^2'  cannot  use 

values  of  © which  would  put  or  02  out  of  the  range  (-it,  +it)  . 

We  see  that  if  ©•  takes  on  values  between  0 and  a.  will 

take  on  values  between  -cp . and  - <p,  . These  are  both  in 

(-H,  +ir)  , so  we  may  consider  any  ©j^,  ©2  as  between  0 and 
(As  a further  note,  if  we  used  ®2  **  *AVG  would 

find  that  some  values  of  ©j^  and  ©2  would  have  to  be  excluded 
because  they  cause  or  02  to  fall  out  of  (-it,  +it).) 

We  can  also  observe  that  although  from  the  point  of  view 
of  determining  the  location  of  point  J it  does  not  matter  if 
we  use  ©j^,  ©j^  + 2it,  etc.,  we  must  remember  that  = © - <p^, 

and  so  adding  a multiple  of  2ii  to  ©^  also  adds  it  to  a^.  We 

have  already  determined  that  values  of  ©^^  between  0 and 
will  cause  to  be  in  (-ir,  +ii)  , so  clearly  we  cannot  add 
multiples  of  2n  to  these  values  of  ©^^  without  putting  the 
out  of  (-tr,  +ir)  . 

Finally,  remembering  the  observation  on  curvature  made 
in  connection  with  Figure  1,  we  find  that  the  total  curvature 
in  Figure  2a  is 


4aj^sinaj^ 


4a2sino2 


+ 


But,  using  the  law  of  sines 


= Dsine2/sin(0j^+92)  and 
D2  = DsinGj^/sin  (0j^+©2)  so  we  have 


4sin(0j^+02)  Oj^sinOj^ 
D ^ sin02 


a,sina_ 

+ ) 
sin0j^ 


Since  0j^  + 02  = snd  = 0^  - cp^,  we  can  now  state  our 

problem  analytically  as  finding  the  minimum  value  of 


4sin$^VG  , t0j^-cpj^]sin[0j^-cpj^] 
D ^ sin02 


l©2“'P2l  sin(02-«p2l 

sln©^  ^ 


subset  to  the  constraint  ®2  “ ^AVG'  ®2 

only  take  on  values  between  0 and  Clearly,  is  a 

special  case,  as  this  implies  0j^  = ©2  - O.  We  will  treat  this 
case  later. 

9i_92 


If  we  define  0 


DIF 


then  we  can  write 
9i  92 

®1  - 9i  - “fAVG  “ ^2  ~ 'Pi  S’  T " ®2  " 'Pi 


'P2  'Pi 

= _£.__±.  _0=_* 

2 2 ®2  DIF 


- 0, 


and  ©2  - 92 


" “^DIF  " ®1- 


(©1-91)  sin(0j^-cpj^)  (02-92)  sin (©2-92) 


sin©. 


sin©. 


^~'^DIF~®2^  ^^D^®1^  ®^"^*DIF"®1^ 


sin0. 


sin0. 


= + cos(-4'jjjp)8in(-e2))/sine2 

+ (sin('»jjjp)cos(-ej^)  + cos(4'pjp)sin(-ej))/8in©j^ 


•"  < *DIF‘®1^ ( sin«DIF°°*'®l'°°®*DIF^ 


= Sin4>jjjp[(0j,jp+e2)cote2  + (♦j^jp-e^)  cote^l  + ♦avg^^^^W 


4 8 in# - 

The  total  curvature  can  thus  be  written  as  ;r X the 


above  expression.  If,  for  the  moment,  we  ignore  the  special 
cases  in  which  “^avG  “^DIFF  zero,  we  see  that  this 

function  goes  to  infinity  as  either  ©j^  or  ©2  approach  zero. 
Since  the  minimum  curvature  obviously  does  not  occur  at  either 
end  of  the  range  in  which  ©j^  and  ©2  must  lie,  we  must  look  for 
a local  minimum.  The  derivative  of  the  above  expression  is 


Since  we  are  assuming  that  and  ♦jjjp  are  nonzero, 

4sin4>^ygSin'I>jjjp/D  is  nonzero,  and  we  must  have 


Since  + ©2  = 


*DIF*»2  ^ -■‘dIF^® 


2 ' 2 
sin  ©2  sin  0 


d02 

AVG'  d©~  “ must  obtain 

— - cot0^  - cot02  = 0 . 

1 


t 

I 


r 


I 

I 

1 


i 


I 


Multiplying  through  by  sin0j^sine2  gives 


sine. 


sine. 


sine. 


sine. 


- cosej^sine2  “ cose28inej^ 


sine. 


sine. 


“dIP*®2' 


sine. 


sine. 


- “"♦avg  * “ ■ 


with  the  derivative  in  this  form,  we  can  see  that  as  e^^  -*•  0, 
the  derivative  goes  to  infinity  with  the  s£une  sign  as  “♦jjjpp 
and  as  02  0,  the  derivative  goes  to  infinity  with  the  seune 

sign  as  (Remember  that  0^  and  02  must  have  the  seune 

sign,  so  sinej^/sine2  and  sine2/sinej^  > 0).  Since  we  know  that 
0^^  and  02  must  be  between  0 and  and  we  know  the  behavior 

of  the  derivative  at  these  points,  it  is  easy  to  solve  the 
above  equation  numerically  by  bisection. 

We  still  have  two  special  cases  to  take  care  of.  If 
'*’diF  ~ then  = 92*  this  case,  the  contour  is  just  a 

single  circular  arc.  This  can,  of  course,  be  viewed  as  two 
circular  arcs  of  equal  radius  joined  at  any  point  we  choose 
along  the  contour.  Note  that  the  derivative  is  zero,  as  ex- 
pected, since  the  contour  and  curvature  remain  the  seune  regard- 
less of  where  we  choose  to  join  the  arcs.  If  is  nonzero, 

we  must  still  have  = 0.  In  this  case  = ~q>2*  Since 

0^1^  and  02  must  lie  between  0 and  0^^  = 02  = 0,  so  the 

arcs  must  join  at  some  point  on  a line  between  and  P2 
(Figure  6) . The  curvature  is 


cp.sincp.  -cp,  sin-cp.  , , 

D—  + 'Pl^i'^'Pl^,  dJ  °1 


L 


and  n2  must  be  between  0 and  D.  It  is  not  difficult  to 
show  that  the  minimum  value  occurs  at  ■ D2  ■ D/2.  (Note 
that  since  ^ 0»  0.)  We  thus  obtain  the  result 

we  would  expect  due  to  symmetry.  We  could  also  arrive  at 
the  last  result  by  examining  the  behavior  of  the  expressions 
we  derived  for  the  curvature  and  its  derivative  as  -*• 
although  this  approach  is  more  complicated. 


3 . Cubic  Interpolation 

We  can  take  another  approach,  using  polynomials.  Suppose 
we  have  two  points  and  ?2»  vectors  and  V2r  as 


in  the  diagremv  below. 


We  can  define  a curve  (x(t),  y(t))  which  passes  through  point 
with  vector  derivative  and  passes  through  point  P2 
with  vector  derivative  V2  in  the  following  way; 


and  let  parameter  t 6 [0,1] 


Then  we  want 


' *^1,'  y'(0)  ’ ''l,  ‘ y(l)  ^ ^2,^  y'(l)  ’ 


We  can  satisfy  the  above  conditions  by  defining 


x(t)  = A^t^  + B^t^  + C^t  + 


y(t)  = Ayt^  + Byt^  + Cyt  t Dy 


where 


"2X  * '’lx  * ''2y 


^"■ix  - ^'■ix  - '’2X 


2P 


ly 


3P 


2y 


ly 


ly 


3P 


ly 


^^27 

''2y 


Of  course,  since  we  are  really  only  interested  in  matching 
the  directions  of  and  V2  at  and  Pj#  the  magnitudes  of 
and  ^2  free  pareuneters. 

Note  that  the  solution  curves  are  invariant  with  respect 
to  rotation  and  translation.  This  follows  from  the  following 
two  properties  of  the  solution  curves  generated  by 

,.(t)  , , 

' y{t)  ' ' A t"*  + B t + C t + D ’ ' 

^ y y y y 


1)  This  class  of  curves  is  obviously  closed  under  ro- 
tation and  translation. 

2)  Any  pair  of  points  and  vectors,  ]^2»  ^1'  ^2 

generates  a unique  solution  curve. 


In  other  words,  suppose  we  have  points  Pj^,  P2  and  vectors 
Vj^,  V2»  and  the  corresponding  solution  curve  is  (x(t)  , y(t)). 
Then  if  we  apply  a rotation  and  translation  to  Pj^,  P2,  ^2 

to  get  P^,  P^,  v|,  clearly  the  curve  (x(t),  y(t))  taken 

through  the  seune  rotation  and  translation  will  yield  a solution 
curve  to  the  problem  with  Pj^,  P^#  v^.  But  since  P|,  P^, 

determine  a unique  solution  curve  (x'(t),  y'(t)),  this 
curve  must  be  a rotated  and  translated  version  of  {x(t) , y(t)). 


In  order  to  implement  this  method,  it  was  of  course 
necessary  to  devise  a means  of  determining  the  two  free  para- 
meters (the  magnitudes  of  the  vectors  and  V2)  . The  method 
used  to  determine  these  parameters  is  based  on  the  following 
observation.  Consider  the  figure  below: 


\ 

I 


L is  the  perpendicular  bisector  of  line  segment  ^ 

is  a circular  arc.  It  is  easily  shown  that  a polynomial 
based  completion,  generated  by  the  technique  described  pre- 
viously, will  pass  through  point  x (the  seune  point  through 
which  a circular  completion  passes)  when 

1^1  I = IV2I  - l + cog© 

The  resulting  polynomial  completion  closely  approximates  a 
circular  completion  until  9 becomes  much  greater  than  7r/2. 
Therefore,  to  facilitate  a comparison  between  the  two  methods 
of  generating  completions,  it  was  decided  to  determine  the 
magnitudes  of  and  V2  in  such  a way  that  in  the  symmetrical 
case  diagreunmed  above  [vj^l  = IV2I  = 2D/(l+coscp)  . There  are 
obviously  many  ways  to  do  this.  The  one  used  in  the  progrcim 


4. 


Examples 

To  generate  examples  on  which  to  run  the  programs,  we 
let  both  cpj^  and  cp2  take  on  values  in  multiples  of  v/A.  This, 
of  course,  results  in  eight  values  for  both  and  if 2,  and 
yields  64  cases.  Many  of  these  cases  are,  however,  reflec- 
tions of  other  cases.  Therefore  only  a set  of  cases  which 
are  unique  up  to  reflection  are  shown.  (Of  course,  the 
case  = TT,  (^2  = n is  not.shown,  as  both  techniques  go  to  in- 
finity when  this  situation  is  approached.) 
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Figure  4b. 
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Figure  4d,  cont'd 
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Figure  5a. 


Figure  5c 
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Polynomial  completions. 
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